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' Abstract. Final state medium-induced gluon radiation in ultradense nuclear matter is shown 

£\J . to favor large angle emission when compared to vacuum bremsstrahlung due to the suppression 

of collinear gluons. In a simple model of momentum transfers from the medium to the jet+gluon 
system the phase space distribution of non-Abelian bremsstrahlung is examined to first order in 
opacity. Perturbative expression for the contribution of the hadronic fragments from medium- 
induced gluons to the back-to-back particle correlations is derived. It is found that in the limit of 
large jet energy loss gluon radiation determines the yield and angular distribution of \A(p\ > 5 
di-hadrons to high transverse momenta pr 2 of the associated particles. Clear transition from 
enhancement to suppression of the away-side di-hadron correlations is established at moderate 
Pt 2 and its experimentally accessible features are predicted versus the trigger particle momentum 
f***^ . pTi • The measured away-side correlation function is shown to be sensitive to the broad rapidity 

""■-■^ ' distribution of di-jets and their vacuum and medium-induced acoplanarity. 

a: 
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1. Introduction 

The discovery of jet quenching - the suppression of large transverse momentum hadron 
production in nuclear collisions relative to the expectation from p+p reactions scaled by the 
number of elementary nucleon-nucleon interactions - is arguably the most exciting new result 
from the Au+Au experimental program at the Relativistic Heavy Ion Collider (RHIC). The 
phenomenon of jet quenching has been established via the attenuation of the single inclusive 
particle spectra j2j and the suppression of the back-to-back di-hadron correlations [HI E]. It has 
been interpreted as critical evidence for large parton energy loss El E] in ultradense quark- 
gluon plasma (QGP) - the deconfined state of matter predicted by quantum chromodynamics 
(QCD). 

So far, the observable effects of the medium-induced gluon radiation itself were considered 
to be modest. For tagged jets, assuming thermalization of the lost energy calculated in |Hj, 
transport models predicted an increase in the multiplicity of associated particles limited to 
Pt 2 — 500 MeV [§]. Using the angular gluon distribution from lljj, jet cone broadening was 
found to be < 10% and challenging to detect experimentally even at the CERN Large Hadron 
Collider (LHC) Q31- 

In these proceedings we follow Ref. |12j and demonstrate that a mechanism, based on the 
destructive interference of color currents from hard and soft parton scattering, can ensure 
broad in angle and frequency final state medium-induced emission spectrum in QCD. For large 
energy loss, perturbative fragmentation of the radiative gluons is found to give a dominant 
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Figure 1. (a) Direct diagrams that contribute to first order in opacity n = L/X. (b) Double 
Born diagrams that interfere with the vacuum bremsstrahlung (c) to the same order in ft. 

contribution to the yield of away-side di-hadrons and to significantly alter their correlations at 
transverse momenta much higher than naively anticipated. We note that the medium-modified 
jet topology can be distorted in two particle back-to-back measurements by the broad di-jet 
rapidity distribution and the vacuum and medium-induced acoplanarity. If present, it should be 
more cleanly detectable on an event-by-event basis and through jet energy flow measurements. 



2. Angular distribution of the medium-induced non-Abelian bremsstrahlung 

We first recall that for hard perturbative scattering the radiative spectrum of real gluon emission 
for small and moderate frequencies u is given by |Hj 
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The distribution in Eq. (^Q) corresponds to the diagram in Fig. ^ where the blob represents the 
large virtual momentum exchange. In Eq. fl* = arcsin(A;_|_/u;) is the angle relative to the 
jet axis, 5 is the azimuthal cone angle (both illustrated in Fig. |2J), a s is the strong coupling 
constant and Cr = 3 (4/3) for gluon (quark) jets, respectively. Virtual gluon corrections remove 
the uj — > infrared singularity in the cross sections in accord with the Kinoshita-Lee-Nauenberg 
theorem ^3] but the collinear 6* — > divergence has to be regulated or subtracted in the parton 
distribution functions (PDFs) and the fragmentation functions (FFs). 

In contrast, the final state medium-induced bremsstrahlung spectrum is both collinear and 
infrared safe. The full solution for the medium induced gluon radiation off jets produced in a 
hard collisions at early times Tj e t — 1 / E inside a nuclear medium of length L can be obtained to 
all orders in the correlations between the multiple scattering centers via the reaction operator 
approach [S] 
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where ^2 = is understood. All vectors in Eq. (0) are two dimensional and coincident with 
the plane orthogonal to the direction of jet propagation. We denote by 6 the unit vector in 



the direction of the radiative gluon's transverse momentum, k = u sin 9*5. The color current 
propagators read 
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The momentum transfers are distributed according to a normalized elastic differential cross 
section 
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which models scattering by soft partons with a thermally generated Debye screening mass 
fJ,(i). From Eq. (J3J) (q|) oc p 2 {zi) and for a quark- gluon plasma in local thermal equilibrium 
/U 2 (zj) ~ 4Tra s T 2 (zi). It has been shown [S] via the cancellation of direct and virtual diagrams 
that in the eikonal limit only the gluon mean free path X g (i) enters the medium-induced 
bremsstrahlung spectrum in Eq. (j2J. For gluon dominated bulk soft matter a e i(i) ~ ^na 2 / /J, 2 (i) 
and \ g (i) = l/a e i(i)p(i). For the case of (1+1)D dynamical Bjorken expansion of the QGP 0E] 

Ufa) = mOo) y—j , A ff (^) = X g (z ) i^j-j . (5) 

The characteristic path length dependence, oc L 2 for static plasmas, of the non-Abelian energy 
loss in Eq. (j2J comes from the interference phases and is differentially controlled by the inverse 
formation times, 

(uj sin 9*5 - q m q n ) 2 

w (m,---,n) = 2^ ) (OJ 

and the separations of the subsequent scattering centers Az^ = z^ — Zk—i- It is the non- 
Abelian analogue of the Landau-Pomeranchuk-Migdal destructive interference effect in QED |14j . 
The extension to the case of massive quarks can be carried out systematically to all orders in 
opacity JH] by identifying the heavy quark mass terms oc uj 2 M 2 /E 2 . 

It can be shown from Eq. (j2J that the propagator poles, Eq. Q, where potential divergences 
may arise are precisely canceled by the interference phases, Eq. ©, to all orders in opacity. 
This point will be better illustrated if we examine the Gyulassy-Levai-Vitev [S] (GLV) gluon 
distribution in angle and frequency to first order in the mean number of soft interactions in the 
plasma. The diagrams shown in Fig. ^ yield: 
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In Eq. (JJJ) the two terms, |-M a | 2 and 2Re7W^7W c are individually divergent. Their sum is, 
however, well defined due to the cancellation of the soft and collinear singularities. From Eq. (J7J) 
the gluon distribution is not only finite when 9* — > but vanishes on average due to the 
uniform angular distribution of momentum transfers from the medium, da cos a = 0. Here 
cos a = qi • <5/|qi|. We have checked that for physical gluons of |k| < uj the cancellations 
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Figure 2. The angular distribution of medium-induced bremsstrahlung of E = 6 GeV gluon 
jet for fixed values of the radiative gluon energy u> = 0.5, 1,2,4 GeV. Top and bottom panels 
represent (1+1)D Bjorken expanding medium of transverse size L = 9 fm and L = 3 fm, 
respectively. Inserts show (8*) versus uj. Right panels illustrate the gluon spectrum and radiation 
intensity of E = 40 GeV gluon jet. 



discussed here persist to all orders in the mean number of scatterings jB]. The small frequency 
and small angle spectral behavior of dJS ^ icd /duod sin 6* d5 remains under perturbative control. 

We also emphasize that destructive quantum interference suppresses radiation of Az <C If (Hj • 
Gluons of final transverse momentum k are initially much more collinear with the jet but acquire 
their moderately large angle 8* from momentum transfers from the medium. The formation 
length at the emission vertex, which allows for interference on the scale of the nuclear size, reads 

l ^ 2a; 

/ (k - qi) 2 q 2 _ 2oo sin 0* qi • 5 + to 2 sin 2 8* ' 

The induced radiation decouples from the jet at It ~ Az oc L/2 and reduces the possibility of 
gluon showering due to further interactions in the plasma. For example, in central L = 9 fm 
collisions (If) ~ 4 f m and with formation times zq = 0.25 — 0.5 fm by the time a physical gluon 
is formed the plasma density has dropped by a factor of 8 to 16 from the longitudinal Bjorken 
expansion. 

Given the vastly different angular behavior of the vacuum and the medium-induced gluon 
bremsstrahlung, Eqs. Q and l[7|). it is critical to identify the phase space where cancellation 
of the color currents induced by the hard and soft scattering occurs. We fix the parameters 
of the medium in Eq. to h(zq) = 1.5 GeV and X g (zo) = 0.75 fm at initial formation time 
zq = 0.25 fm. Since small |k| and u> emission is suppressed, we use only a moderate a s = 0.25. 
Triggering on high pr x hadron directs its parent parton "c" away from the medium and places 
the collision point of the lowest order (LO) ab — ► cd underlying perturbative process JH] close 
to the periphery of the nuclear overlap region. Then, it is the back-scattered jet "d" that 
traverses the QGP. For large nuclei, such as Au and Pb, path lengths L = 9 (3) fm are used to 



illustrate central (peripheral) collisions, respectively. We limit gluon emission to the forward jet 
hemisphere, < 9* < %. 

The angular distribution of medium-induced radiation for E = 6 GeV gluon jet for select 
values of u is shown in the right panels of Fig. [2j We find that gluon emission is strongly 
suppressed within a cone of opening angle 9* ~ 0.25 rad due to the cancellation of collinear 
bremsstrahlung - a mechanism different from a Gaussian random walk in 9*. The broad gluon 
distribution can be characterized by the mean emission angle 

-l 

(9) 

given in the inserts of Fig. [2j Left panels illustrate the infrared safety property of the spectrum. 
Energy loss is dominated by semihard few GeV gluons, see ujdN 9 / doj. The gluon number 
dN 9 /duj distribution peaks at even smaller frequencies uj. Momentum transfers from the medium 
|qi| ~ 1 — 2 GeV are clearly able to generate broad angular bremsstrahlung distribution. 
However, it is the cancellation of the 9* — > radiation together with the model of the medium 
which ensure the weak uj dependence of the mean emission angle {9*) ~ 0.5 rad seen in Fig. [3 
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3. Suppression of the single inclusive hadrons 

Dynamical nuclear effects in multi-particle production in A+B collisions can be studied through 
the ratio 

R (n) = da%' hn / d Vl ■ ■ ■ dy n d 2 p Tl • • • d 2 p Tn 
AB (N c A °g) da h N ^ h "/d yi ■ ■ ■ dy n d 2 PTl ■ ■ ■ d?p Tn ■ 

In Eq. (0 {N$j) is the mean number of binary collisions which can be calculated from an 
optical Glauber model. 

Results form the perturbative calculation of the nuclear modification to the neutral and 
charged pion production in central Au + Au collisions are given in the left hand side of Figure El 
At all energies there is cancellation between the Cronin enhancement, which arises from the 
transverse momentum diffusion of fast partons in cold nuclear matter j5J 1171 1181 119j , and the 
subsequent inelastic jet attenuation in the final state. This interplay is most pronounced at 
the low Super Proton Synchrotron (SPS) s/snn = 17 GeV where, in the absence of quenching, 
the corresponding enhancement could reach a factor of 3-4 With final state energy loss 
taken into account, R^ A (p Tl ) is shown versus the WA98 [201 and the CERES [21] Pb+Pb and 
Pb+Au data with p+p baseline defined as in [221- Systematic comparison of the calculated 
nuclear modification to the existing experimental measurements can also be found in [22] • The 
agreement between data and theory is good at high pr 1 but the quenching is over predicted at 
PTx ~ few GeV. 

In two particle correlation measurements the energy loss of the second parent parton at high 
Pt 17 Pt 2 qualitatively leads to 1 < ^aaI^aa — Numerical estimates in Fig. Etb) from jet 
quenching alone indicate that the double inclusive hadron suppression is 25% — 40% larger than 
the single inclusive quenching in the 5 < pt x = Pt 2 < 7 GeV range at y / s/v]v = 62 GeV. Double 
inclusive cross section quenching is manifest in the attenuation of the away-side correlation 
function C 2 (A0) = (l/N trig )dN hlh2 /dAcj) 23 j. In Fig^c) it leads to a 4 - 5 fold attenuation of 
the area ^4f ar relative to the p+p case. In contrast, elastic transverse momentum diffusion will 
only result in a broader C2(A(/>). We note that the experimentally measured suppression value 
will be sensitive to the subtraction of the elliptic flow V2 component which may remove part or 
all of the residual correlations. 

To study the details of non-Abelian energy loss in the QGP we extend the calculation of the 
quenching of the away-side di-hadrons to low values of pt 2 hi the next section. 
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Figure 3. (a) Calculated nuclear modification factor Raa{pt x ) versus the center of mass 
energy y / ijv r /v = 17, 62, 200 GeV for central Au+Au collisions. The yellow band represents 
the perturbative QCD expectation for the depletion of the neutral pion multiplicity at ^/s^n = 
62 GeV. (b) The ratio of the attenuation of the single and double inclusive tt° cross sections 
Raa(PTi) / RjIa 2 (PTi = Pt 2 ) m central Au+Au collisions at ^/snn = 62 GeV and (c) the 
manifestation of the double inclusive hadron suppression in the quenching of the away-side 
di-hadron correlation function C2(A<^>) 23 . 



4. Redistribution of the lost energy 

Results reported in Ref. ^2] are important for the future LHC heavy ion program since they for 
the first time suggest a large, and therefore detectable, broadening of the jet cone in the nuclear 
environment. Note that {9*) implies a redistribution of the energy flow even within jets of large 
opening angle R = y At] 2 + A(j) 2 = 1. If high statistics measurements of jet quenching away 
from midrapidity become available at RHIC II novel "octupol twist" moments of the parton 
absorption pattern in the QGP may be accessible [2Hj. At present, however, a key question for 
perturbative QCD phenomenology is whether the medium induced gluon bremsstrahlung can 
significantly alter the di-hadron correlations measured at midrapidity at RHIC |H1 IU 1241 125j . 
We naturally focus on the away-side \A(p\ > ^ case, where medium effects are the largest. 
Nuclear modifications build upon the LO double inclusive hadron production cross section, 
which is calculable in the perturbative QCD factorization approach [23| if either of the hadrons 
is moderately hard (pn or px 2 > few GeV) 123] 

da NN is f 1 dz l n / \ [ 7-) / WA n! <t>a/N(Xa)<l>b/N{xb) 
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x 2ira 2 s \M ab ^ cd \ 2 . (11) 

In Eq. (fTTj) K = 2 is a next-to-leading order K-factor, x a ,b = Pa,b/PN a ,N b are the momentum 
fractions of the incoming partons and Z\% = Ph 1 ,h 2 /Pc,d are the momentum fractions of the 



hadronic fragments. We use standard lowest order Gluck-Reya-Vogt PDFs [2B] and Binnewies- 
Kniehl-Kramer FFs jSHj- Renormalization, factorization and fragmentation scales are suppressed 
everywhere for clarity. The spin (polarization) and color averaged matrix elements — >cd\ 
are given in |16j . 

We first recall the physical effects that alter the LO perturbative formula, Eq. IjllJ) . A 
modification that does not change the Ac^-integrated cross section is vacuum- and medium- 
induced acoplanaritv jlHl I24j . The deviation of jets from being back-to-back in a plane 
perpendicular to the collision axis arises from the soft gluon radiation and transverse momentum 
diffusion in dense nuclear matter |18j . In the approximation of collinear fragmentation, the 
width of the away-side hadron-hadron correlation function can be related to the accumulated 
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a good description of \A(p\ > \ correlations measured in elementary p+p collisions at v S 
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quark (and a 2.25 larger value for a scattered gluon). Additional broadening arises from the 
interactions of the jet in the QGP that ultimately lead to the reported energy loss. Using the 
parameters of the medium for (1+1)D expansion we find 



(^T hot) 
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although only half is projected on the </?-plane, (/c^)^ = (kj, vac ) v + ^(kj^ hot ). 

Two competing mechanisms do, however, change the px 2 dependence of the perturbative cross 
section, Eq. (|11|) . First is the the parent jet "d" fractional energy loss e = AE^/Ed, which we 
here for simplicity consider on average and evaluate by integrating Eq. Q • It leads to a rescaling 
of the hadronic fragmentation momentum fraction Z2 — > ^2/(1 — e) Physically, fewer high 
Pt 2 particles are produced by the attenuated parton of energy Ed — AEd- If the energy loss 
is large, a second mechanism is invoked as a consequence. Hadronic fragments of the radiative 
gluons will increase the probability of finding low and moderate pt 2 particles associated with 
the interacting jet [2]. 

To calculate di-hadron correlations, we first map the jet structure of a hard 90°- scattered 
parton on rapidity y ~ 77 = — lntan(#/2) (9 being the angle relative to the beam axis) and 
azimuth <fi, tan 2 9* = cot 2 9 + tan 2 <j> , tan 5 = — cot#/tan(/>. The approximately flat rapidity 
distribution of the away-side jet near 2/2 = can be used to sum over all emission angles 
e (6> mi n, #max) C (0, 7r) yielding 
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The Jacobian of the transformation in Eq. (H 
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It is critical to note that projection on a plane coincident with the jet cone axis (the 0-plane in 
Eq. (|14[) is one example) efficiently masks the 9* — > void in the angular distribution of medium- 
induced gluons reported in Fig.|2j While projecting the medium induced radiation on the y = 




Figure 4. Nuclear modification of the back-to-back di-hadron correlations with and without 
the contribution of medium-induced bremsstrahlung. Top and bottom panels illustrate central 

(2) 

and peripheral collisions, respectively. Data is from STAR 4 j. Right panels predict R\ A for low 
and high pji hadron triggers. 



plane is an approximation with the largest deviations being close to the edges of phase space 
(in this case ^ and t£ ), it is much less expensive computationally than a full simulation of the 
away-side jet+gluons. Our conclusion about the importance of the broad rapidity distribution 
of di-jets and the vacuum and medium-induced acoplanarity is independent of the physical 
mechanism that depletes the parton (or particle) multiplicity in a cone around the jet axis. 

The end analytic result for the modification to Eq. ()11() per average nucleon-nucleon collision 
in the heavy ion environment can be derived from the energy sum rule for all hadronic fragments 
from the jet, 
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Here, z g = Pt 2 / uj - 

Whether medium-induced gluon radiation may have significant observable consequences for 
the large angle di-hadron correlations depends on its relative contribution to the Ac/?-integrated 
cross section. From Eq. 1)11(1 this can be studied versus pt 2 via the ratio RJa (23j . Numerical 
results, shown in Fig. 21 correspond to triggering on a high px 1 =4-6 GeV pion and measuring 
all associated 7r + + 7r° + 7r~. Depletion of hadrons from the quenched parent parton alone leads 
to a large suppression of the double inclusive cross section with weak px 2 dependence. Hadronic 
feedback from the medium-induced gluon radiation, however, completely changes the nuclear 

(2) 

modification factor R\^. It now shows a clear transition from a quenching of the away-side 
jet at high transverse momenta to enhancement at pt 2 < 2 GeV, a scale significantly larger 
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Figure 5. The angular distribution of |A</?| > ^ di-hadrons for central and peripheral A+A 
collisions. Solid and dashed lines give separately the contribution of the radiative gluons and 
the attenuated parent parton. The lower right panel compares the the width of the correlation 
function in p+p to the one of the suppressed parent parton in central A+A reactions. 



than the one found in j^j. In fact, we have checked that for large AEd the back-to-back di- 
hadron correlations are dominated by radiative gluons to unexpectedly high px 2 ~ 10 GeV. 
It should be noted that the experimental trigger introduces a bias toward large z±, which is 
difficult to account for analytically. Therefore we expect deviations between theory and the 
measurement when px 2 > Pt±- STAR data in central and peripheral Au+Au collisions [1] is 

(2) 

shown for comparison. Experimentally testable predictions for the shape and magnitude of R\ A 

(I} ■ 

and the enhancement-to-suppression transition R AA = 1 versus p r r 1 are also given in Fig. 0] It 
has been previously argued [30] that detailed balance reduces jet suppression at partonic scales 
< few//D, though such mechanism cannot produce enhancement of the low transverse 
momentum particle production observed in Fig. |I] 

With a dominant contribution to the di-hadron yields, the medium induced gluons are bound 
to also determine the width of the correlation function. Two-particle distributions in A+A 
reactions, calculated from Eqs. and (jlfij) for a pj\ = 6 GeV trigger pion and two different 
Pt 2 = 1, 4 GeV associated pions, are shown in Fig. ©. Qualitatively, the medium-induced 
gluon component to the cross section controls the growth of the correlation width in central 
and semi-central nuclear collisions. Quantitatively, the effect should be even larger than the 
one estimated here, which is limited by the imposed < 8* < \ constraint. The decisive role 
of the bremsstrahlung spectrum, Eq. 0, in establishing the Ac^-shape of away-side di-hadrons 
is further clarified in the bottom right panel of Fig. [3 Hadronic fragments from the quenched 
parton "d" only are shown to yield a distribution that is not broader than the one anticipated 
in p+p reactions. Experimental measurements of significantly enhanced widths for |A<^| > l| 
two-particle correlation in A+A collisions should thus point to copious hadron production from 
medium-induced large angle gluon emission. 



5. Conclusions 

In summary, we calculated the transverse momentum and angular distribution of the away- 
side di-hadron correlations in the framework of the perturbative QCD factorization approach, 
augmented by inelastic jet interactions in the quark- gluon plasma. At RHIC energies we found 
that the medium-induced gluon radiation determines the two-particle yields and the width 
of their correlation function to surprisingly high transverse momentum px 2 ~ 10 GeV. Clear 
transition from back-to-back jet enhancement to back-to-back jet quenching is established at 
moderate pt 2 = 1 — 3 GeV, independent of collision centrality but sensitive to the trigger 
hadron momentum p^ x . Definitive experimental determination of its features, predicted in 
Ref . , will for the first time provide a handle on the differential spectrum of medium- induced 
non-Abelian bremsstrahlung. Coincidental confirmation of large broadening of the away-side di- 
hadrons would require a critical reassessment of the origin of intermediate transverse momentum 
particles in central and semi-central nuclear collisions. A full numerical simulation of jets away 
from midrapidity is required for more quantitative studies of the two particle topology. For jet 
physics at the LHC, our findings ensure a measurable increase in the jet width in ultradense 
nuclear matter. 
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